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^ Abstract 

We study large scale structure in the cosmology of Coleman-de Luccia bubble collisions. 
Within a set of controlled approximations we calculate the effects on galaxy motion 
O seen from inside a bubble which has undergone such a collision. We find that generically 

>■ bubble collisions lead to a coherent bulk flow of galaxies on some part of our sky, the 

details of which depend on the initial conditions of the collision and redshift to the 
;JJ galaxy in question. With other parameters held fixed the effects weaken as the amount 

^ of inflation inside our bubble grows, but can produce measurable flows past the number 

of efolds required to solve the flatness and horizon problems. 
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1 Introduction 



One of the greatest challenges for string theory is to find a testable prediction which can 
differentiate it from other potential theories of quantum gravity and high energy physics. 
Unfortunately, the typical energy scales at which string theoretic effects become important is 
likely far out of reach for any earth-based particle accelerator. Cosmology though, provides 
another opportunity. The early universe was arbitrarily hot and dense and signatures of the 
physics at that time remain in the large-scale structure of the universe today. 

It has become increasingly clear in recent years that one of the characteristic features 
of string theory is the "landscape" model of cosmology [1, 2, 3]. There are many ways to 
generate realistic four-dimensional universes in string theory via flux compactifications. Var- 
ious parameters specifying the compactification become scalar fields in our four-dimensional 
universe. The landscape is a potential energy function of the would-be moduli scalar fields, 
which possesses many distinct minima. These minima are meta-stable and can decay via 
bubble nucleation-type tunneling or quantum thermal processes. In the early universe, one 
expects all the minima to be populated as the temperature decreases and the fields begin 
to settle. The minima with the highest values of potential energy then begin to inflate, 
and are expected to rapidly dominate the volume. After a few string times the universe 
will be almost completely full of rapidly inflating regions, and as time continues to pass the 
metastable field values in these regions will begin to decay by bubble nucleation. Since we 
live in a region with very small vacuum energy, we expect that our observable universe is 
contained in such a bubble. 

The most promising way to test this idea is to look at large scale cosmological structure. 
If we model the string theory landscape as a scalar field coupled to gravity with a potential 
with several minima, universes inside a bubble which formed via a tunneling instanton have 
certain characteristic features: they are open, they are curvature dominated at the "big bang" 
(rendering it non-singular), and they have a characteristic spectrum of density perturbations 
generated during the infiationary period that takes place inside them [4, 5, 6, 7, 8, 9, 10]. 

Moreover, because the bubble forms inside a metastable region, other bubbles will appear 
in the space around it. If these bubbles appear close enough to ours, they will collide with 
it [11, 12, 13]. The two bubbles then "stick" together, separated by a domain wall, and 
the force of the coUision releases a pulse of energy which propagates through each of the 
pair. This process (as seen from inside our bubble) makes the universe anisotropic, affects 
the cosmic microwave background and the formation of structure. A collision between two 
bubbles produces a spacetime which can be solved for exactly [14, 15, 16, 17]. 

In [15, 17], one of these authors began the study of this process and its effects. In 
these papers we principally studied the effects of such a collision on the cosmic microwave 
background (CMB). We found that bubble collisions naturally create hemispherical power 
asymmetries spread through the lower CMB multipoles and/or produce relatively small cold 
or hot spots, depending on the time of the collision (and hence the size of the disk), and 
can lead to measurable non-Gaussianities. The qualitative features, however, are always 
the same: a collision with a single other bubble produces anisotropics that depend only on 
one angle — the angle to the vector pointing from our location to the center of the other 



2 



bubble — and which are generically monotonic in the angle. 

While the CMB temperature map is a good probe of these effects, it is by no means the 
only one. In this paper, we begin to look for other measurable effects by studying how a 
single bubble collision affects large scale structure and galaxy motion. We find that collisions 
generically lead to a large scale coherent flow of galaxies in our sky. The size of the affected 
disk of galaxies, and the magnitude of the effect are controlled by similar order parameters 
for the CMB effect, as well as the redshifts of the galaxies in question. Generically, the 
angular size of the affected disk on the sky for large scale structure is always larger than 
the corresponding disk for the CMB hot/cold spots. It is very interesting to consider our 
results in light of recent observational measurements of bulk galaxy flows dubbed "dark flow" 
[18, 19, 20]. 

Our model consists of a collision between our bubble - which we treat as a thin-wall 
Coleman-de Luccia bubble which undergoes a period of inflation after it formed (which we 
model as a de Sitter phase) followed by a period of radiation domination - and another 
different bubble, the details of which are largely unimportant. We assume that the tension 
of the domain wall between the bubbles and the vacuum energy in the other bubble are 
such that the domain wall accelerates away from us. Such collisions can produce observable 
effects which are not in conflict with current data. We will also assume that the inflaton 
is the field that underwent the initial tunneling transition, so that the reheating surface is 
determined by the evolution of that field and can be affected by the collision with the other 
bubble through its non-linear equation of motion, and that no other scalar fields are relevant. 

After some number of efolds of inflation, our bubble will reheat and become radiation 
dominated. The subsequent expansion and the geometry of the reheating surface itself will 
be modifled by the presence of the domain wall. In [17], it was shown how the surface was 
modified. This modified surface will also backreact on the spacetime because the various 
fiuids created at reheating will no longer be classically comoving. In [17] this backreaction 
was ignored and it was assumed the universe could be approximated as standard radiation 
dominated Fricdmann-Robertson- Walker (FRW) universe for the subsequent evolution. We 
will show here that this approximation is valid for the CMB photons, but not for the matter 
fiuid that will form the galaxies we are interested in here. To find the evolution of galax- 
ies, some measure of this backreaction must be taken into account, and we make a first 
approximate attempt to do so here. 

There are a number of anomalies in current cosmological data, including a cold spot 

[21, 22, 23, 24, 25], hemispherical power asymmetry [26], non-Gaussianitics (see e.g. [27]), 
the "axis of evil" [28], and the tentative possibility of a "dark flow" [18, 19, 20]. While the 
signiflcance of these effects is unclear, there are few if any flrst-principle models which can 
account for them. Our results here, combined with [15, 17] indicate that bubble collisions can 
naturally create hemispherical power asymmetries spread through the lower CMB multipoles, 
relatively small cold or hot spots, measurable non-Gaussianities and a possible macroscopic, 
coherent bulk flow of galaxies. While it is not yet clear that a single collision can produce 
all of these effects in agreement with observations, it is nevertheless interesting that bubble 
collisions can at least qualitatively produce many of these effects. 
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2 Basic Setup 



In this section we will set up the basic scenario for the bubble collision and subsequent 
evolution we wish to study. The details of the dynamics of general collisions were studied in 
[15], and some effects on the CMB were studied in [17]. We refer the reader there for a more 
complete treatment. We will focus on a collision scenario where our bubble is dominated 
by positive vacuum energy at formation, so that it undergoes a period of inflation. We will 
further assume the pressures are such that domain wall that forms between the two bubbles 
is accelerating away from the center of our bubble at late times. This second assumption 
is valid for collisions between two dS bubbles if ours has a smaller cosmological constant.^ 
While more general scenarios are possible, these assumptions lead to an evolution that is 
consistent with current observations. 

We study the motion of galaxies and large scale structure in this scenario. A full treatment 
would require solving the full cosmological perturbation theory and evolution taking into 
account the macroscopic anisotropy caused by the collision. This treatment is outside the 
realm of the current paper. Instead, we will make several assumptions so that an approximate 
analytic treatment is possible. Our basic scenario unfolds as follows: our bubble undergoes 
a period of inflation sufficient to solve the horizon, flatness and relic problems in line with 
current observations {N > 65 e-folds) after which the spacetime reheats. After reheating 
the universe can be approximately described by a radiation dominated (RD) Friedmann- 
Robertson- Walker (FRW) universe. We say approximately, because as shown in [17] the 
reheating surface is modified by the bubble collision and formation of the domain wall. In 
short, the domain wall causes the reheating surface to "bend" towards or away from the 
domain wall depending on the details of the coUision. See Fig. 1 and [17] for more details. 
Because the universe now reheats along a deformed, non-comoving surface the spacetime is 
not truly RD-FRW, as some of the symmetry of this geometry is broken. In [17] we still 
approximated the subsequent evolution as taking place in RD-FRW. For CMB photons, this 
is a reasonable approximation (as will be shown later). However, for the matter fluid and 
galaxies, this approximation is no longer justifled and some measure of the backreaction and 
modified geometry must be taken into account. We will attempt to do so here. In reality, 
the spacetime will also undergo a transition to matter domination at some later time after 
last scattering. We ignore this latter transition in this paper, as we believe the essence of 
the galaxies motion is captured in our setup. It would nevertheless be interesting to study 
this in more detail and see if our assumptions hold. 

To find the approximate motion of galaxies, we assume that the initial conditions of the 
matter fluid are that it is moving orthogonal to the reheating surface. We will then attempt 

to both capture the backreaction of this fluid on the spacetime (thereby modifying our 
spacetime from RD-FRW) and track the motion of galaxies through this new spacetime by 
computing the timelike geodesies for each galaxy. We are thus making several assumptions, 
some of which we outline here: 

^It also occurs when our bubble has a larger cosmological constant or in collisions with AdS bubbles when 
certain conditions involving the tension of the wall are satisfied [14, 15, 16, 29]. 
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• We compute the backreaction by assuming that the local geometry around any point 
on the reheating surface is RD-FRW. However, each patch will not necessarily be 
comoving with respect to other patches and so a global description as RD-FRW is 
impossible. 

• We do not account for further evolution of the backreaction effect, i.e. self interaction 
and gravitation of the radiation fluid and its subsequent effect on the geometry. We 
expect that these will have a negligible effect on the galaxy motion we are interested 
in here. 

• We ignore self-interaction of the matter ffuid/galaxies and instead compute freestream- 
ing geodesies in the modified geometry. 

2.1 Review of geometries 

In this section we will give a brief review of the geometries and establish our notation 
conventions, for more details as well as a full description of the causal structure the reader 
is directed to [15, 17]. 

A single Coleman-de Luccia bubble in four dimensions has an S0{3, 1) invariance, inher- 
ited from the 5*0(4) of the Euclidean instanton [4]. When another bubble collides with ours 
a preferred direction is picked out, breaking the symmetry down to S0{2, 1). As a result one 
can choose coordinates in such a way that the metric takes the form of a warped product 
between two-dimensional hyperboloids {H2) and a 2D space. 

Under the assumption that all the energy in the spacetime remains confined to thin 
shells or is in the form of vacuum energy, one can find the solution describing the collision 
by patching together general solutions to Einstein's equations in vacuum plus cosmological 
constant with this symmetry [15, 14, 30, 31]. The overall result is that the spacetime inside 
our bubble is divided into two regions-one that is outside the future lightcone of the collision 
and can be described using an S0{3, 1) invariant metric (a standard open FRW cosmology), 
and one which is inside, is affected by backreaction, and is most conveniently described 
in the general case using 2D hyperboloids. However, we will find that at the late times 
important in our analysis, it becomes possible to once again use 3D hyperboloids and in our 
approximation 3D flat space (that is the geometry will be approximately locally flat FRW). 

de Sitter space: The general SO {2, 1) invariant spacetime with a positive cosmological 
constant is given by [15] 



+ g{t)dx'^ + t^dHl, 



(2.1) 



git) 



where 



9it) = l + 



e t ' 



(2.2) 
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Radiation Shock 




Bubble Walls 



Figure 1: Sketch of the causal structure of the spacetime in the bubble collision scenario. 
We define tc and regions I, II and III of our bubble in Sec. 3 

and 

dH^ = ^ sinh^ pdip'^ (2.3) 

is the metric on the unit H2 and we identify a; ~ a; + iri. We are most interested in observers 
far from the domain wall, as these see cosmologies most like ours. In that case it was shown in 
[17] that for the collision to be within our past lightcone we require to <^t cosmologically 
relevant times. Therefore, throughout the rest of the paper we approximate ~ giving 
g{t) ~ 1 +t^/£^.^ This is a metric for pure dS space (albeit in unconventional coordinates) 
with radius of curvature i. Hereafter, we set £ = 1; it can easily be restored with dimensional 
analysis. These coordinates, which we denote the H2 coordinates, are related to the SO{3, 1) 
or invariant ones by 

coshr = vT+~t^cos(x — Xq), 
sinhr sinh^ cos^ = vT+~t^sin(a; — Xq), (2.4) 

sinhr cosh^ = tcoshp, 
sinh T sinh ^ sin 9 — t sinh p, 

ip = ip, (2.5) 



^It is important to note that to/t may not be small at early times when some fluctuations left the horizon. 
It would be interesting to investigate what effect this has on the CMB and other observables. 
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where the free parameter xq is used to fix the location of the origin of H3. In these coordinates 
the metric becomes 

dsls = -dr^ + sinhVd//|, with (2.6) 

It was also shown in [17] that an observer can only see an exponentially small part of the 
reheating surface in ^, that is the an observer has access to is very small (we will make 
this more precise shortly). Thus, if we choose xq appropriately we have ^ <^ 1 for any 
observer. Throughout we make this approximation and we thus find 

dHl = di^ + sinh^ i {de^ + sin^ 6 dip^) d^^ + {dO^ + sin^ 6 dip^) = dsl^. (2.7) 

Radiation domination: We will also make use of the RD-FRW metric given by 

dsl^^ = -dr^ + a^{T) dHl ^ -dr'^ + a^{T) ds^^, (2.8) 

a(r) = ^Jc{t+\), (2.9) 

with C a constant.^ 

At reheating we need to match both the scale factor a and Hubble constant H of the two 
metrics [17]. This gives 

T^"" = ^(tanhrf-l), (2.10) 

C = 2coshT;'^ sinhr;*^ = 2a^, (2.11) 

where a,. = a{T^^) is the scale factor at reheating. Note that (2.10) is a gluing relation 
joining the two spacetimes; it should not be interpreted as a coordinate transformation. We 
also note that we are able to use this solution only because to/t <^ 1 for the cosmologically 
relevant times. If this were not the case, the solution would be more complicated. Due to 
the reduced symmetry, the general solution for RD with SO {2, 1) invariance is not known. ^ 



Pcirameters and expansions: It will be useful here to define our various parameters and 
orders of expansion. The number of cfolds of infiation is N = t^^. After reheating the 
temperature redshifts as 1/a so we can define 

e^* = T,/T° = ajar, (2.12) 

where the n subscript refers to the observer's time and this defines N^. In [17] it was shown 
that an observer can see a part of the reheating surface approximately fa 4e^*~^ in size. 

^We find it convenient to brealc slightly with the conventions of [17] by translating trd by ^ above; this 
will simplify subsequent expressions. 

^The problem is roughly equivalent to finding the smooth metric describing a Schwazschild black hole 
embedded in an asymptotically FRW spacetime. 
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As ~ 1 is one curvature radius of the H3, is roughly the number of efolds required 
to solve the flatness and horizon problems. To agree with current observations wc require 
N,N^: ^ 1 and N^, < N. This justifies the previous assumption that ^ 1. In our 
computations we frequently drop subleading terms that are suppressed by e^^ or e~^*, and 
expand all quantities as power series in 

e = e^*-^ < 1. (2.13) 

We will work to first subleading order in e throughout this paper. 

3 Reheating and matter flow 

In this section we review and discuss the geometry of the reheating surface, and its effect on 
the post-reheating spacetime. A full treatment of the reheating surface is given in [17] and 
for details the reader is directed there. Here we will present a brief review. The basic idea is 
as discussed in Sec. 2: the collision and domain wall modify the reheating surface and break 
the 5*0(3, 1) isotropy to 5*0(2, 1). Wc define the reheating surface as a surface of constant 
energy density for the infiaton field (p. During infiation these surfaces are also surfaces of 
constant value for the inflation fleld. This agrees with the usual slow-roll conditions so long 
as we are far enough away from the domain wall. Near the wall, the infiaton develops large 
spatial gradients and the usual slow-roll treatment is not enough. In this paper, we will 
focus on the region where 4> = constant is enough to define the reheating surface.^ 

3.1 The scalar field and reheating surface 

In [17] an approximate solution for the scalar field and reheating surface was presented.^ 
This computation involves several steps. We need to solve the scalar field equation with a 
linear potential, 1^(0) = which is a good approximation during slow roll inflation. 

The remaining inputs are the boundary conditions. We expect the scalar field to have a 
set of values (f)^ in the other bubble, some other set of values 0/? in our bubble and a third 
set in the mctastable vacuum outside both bubbles 0m- In the thin wall limit the scalar field 
jumps from near some other local minimum (either the false vacuum our bubble tunneled 
from, or the vacuum inside the bubble it collided with) to a location a bit up the "hill" 
from our local minimum. This transition occurs all along the domain and bubble walls, so 
the scalar field is constant (though not necessarily the same constant) along each of these 
surfaces. As in [17] we choose without loss of generality = along the bubble wall and 
4> = k along the domain wall. Note that /c = is a possible boundary condition. The 
remaining condition on is that it is continuous across the radiation shell of the collision. 
At linear order this is enough to give us the solution for and the reheating surface as a 
function of the potential and the details of the collision (such as collision time, acceleration 

'''Work on the strongly coupled region is in progress [32] . 

full treatment requires solving the non-linear coupled scalar field and gravity system for the collision 
spacetime. Some work in a special case was presented in [33]. 
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of the domain wall, etc.). Once again for a full treatment the reader is directed to [17], where 
the full solution for x) is given. 

For reasons stated above, we are interested in the reheating surface near the radiation 
shock, as illustrated in Fig. 1. The null line for the radiation shock was found in [17] as 

71 

Xrad = arctant — 2 arctantc ~ 2 ~ ^ arctantc, (3.1) 

where tc is the time of the collision, and we approximated arctan t ^ |, as reheating happens 
after t ~ ^ 1. Expanding the solution for the infiaton 0(t, x) around x — Xrad we find 

(l){t,x) ^ -^ln[tcos(a;-a;rad)] - - x,ad)Q{-{x - x,^d)) + 0{{x - x.^df), (3.2) 

where 0ref = — | In [i cos(a; — Xrad)] corresponds to the solution in the unperturbed case, 
without a collision. The step function ensures the solution is given by ^^ef in the part of 
spacetime causally disconnected from the collision. The coefficient of the linear term is given 
by 

4tc 1 



A=--\tr.---^ 



1+t'i TT— 2 arctan tc 



2 1 tr TT — 2 arctan tr 



(3.3) 



up to corrections of order \ ~ . For reahstic values of the parameters (//, k, tc), the slope 
A is of the order 0{1). In the rest of this paper we will be working to first order in the 
parameter A{x - Xrad) ~ 0{e^*~^) < 1? 

It is useful to convert this to coordinates for use in the rest of the paper. Using the 
coordinate transformation (2.4) with the origin of fixed to the radiation shock by setting 
= Xrad, we find 

0(r,e,^)«^0ref(r)-e(-cos^) [|Aecos^ + 0(a] ln2, (3.4) 

where </)rcf(T) = —^r/di and we can drop the In 2 term without loss of generality as it amounts 
to a constant shift of and has no effect on the reheating surface. Hence in these coordinates 
the reheating surface consists of two fiat 3D-surfaces meeting at an angle, as indicated in 
Fig. 2. 

The form of the scalar field in (3.4) suggests that it will be convenient to switch to a 
cartesian system 

i = ^x"^ + ^ ra: = {cos^, 

Q — arctan |, [ y = ^sin6', 

so that 0, and thus the natural equal time surfaces, only depend on r and x? 



''The coefficient of the quadratic term in (3.2) is also of order and hence the linear approximation 

is justified. 

^This X should not be confused with the x in Hi coordinates; from now on we'll only use H-^ coordinates 
and X will be as defined above. 
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3.2 Matter flow and causal regions 



The spacetime after reheating is separated into three causal regions, denoted I, II and III in 
Figs. 1 and 2. Observers in region I are causally disconnected from the bubble collision, and 
the spacetime in region I is RD-FRW space with time coordinate t/ = r, as given in (2.8). 
Since region II is casually disconnected from region I, observers in region II only see the 
perturbed part of the reheating surface. The metric in region II is sourced by the perturbed 
part of reheating surface only. The symmetry of the reheating surface then implies that 
region II is also given by RD-FRW space, but with a time coordinate tu that is constant on 
this part of the surface. One way to view this region is as a RD-FRW space that has been 
boosted to be comoving with respect to the sloped reheating surface. 

Our first goal will be to seek a relation between the II and I coordinates to find a global 
description of these regions (in fact, we will find these relations closely resemble a Lorentz 
boost). In region III, the metric is unknown. A full treatment would backreact the interacting 
radiation fluid in this region, which is beyond the scope of the current paper. Instead, we 
will attempt to approximate the backreacted metric in region III by taking a well motivated 
interpolation between regions I and II. More details on this interpolation will be given in 
Sec. 3.3.1 

We take the matter fluid to be created at reheating, and moving orthogonal to the 
reheating surface. Hence, in regions I and II the matter fluid moves along geodesies defined 
by drj and drjj respectively. This assumption is justified as any peculiar velocity given to the 
matter field at reheating will quickly redshift away in RD-FRW space, leaving the motion 
along and dr^j in regions I and II to a good approximation. 

3.3 Relations between the regions and metrics 

Region II: From the form of the scalar field (3.4), and the gluing relation (2.10), we find 
that the reheating surface is given in RD coordinates as 

Reheating surface: [^^^^^^ ^ ^ ^""^ (3-6) 

To define the coordinates in region II we note that r/j = on the sloped part of the 
reheating surface (3.6), which together with the space being RD-FRW implies 

Tii = 7(t7 + Axi) (3.7) 

for some constant 7. In a coordinate basis we have (dr//, dxjj) = 0, where (, ) denotes the 
inner product with respect to the metric, and thus 

= ^{Adr, - (3.8) 

for some 7. We can fix the constants of proportionality at the intersection of the two parts 
of the reheating surface, as the metric there is given by RD-FRW in both coordinates, and 
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Figure 2: Spacetimes after reheating, showing the reheating surface, an equal time surface tq, 
and the three time axes. The tilting of the time axes is exaggerated to distinguish between 
them; in reality they would be almost parallel for large a^. The figure on the left (right) 
corresponds to A > (< 0) and would lead to a hot (cold) spot in the CMB. 

we have 

-1 = {drn, drii) = j'^{dTi + Adxj, dri + Adxi) = 7^(-l H — ^), 
Or = {9x11 > 9x11 ) = 7 (^^r/ - d^i , Adrj - ) = 7^ ( -^^ + ) , 

T = T=^=. (3.9) 



This gives the coordinate relations 

Tn = i{ti + Axi) , 

xn = ~^{^^^^^^' {^AQi) 

yii = yi: <Pii = <Pi- (3.11) 
We can write the metric in region II in I-coordinates as 

dsn = -drfj + a\dxjj + dyjj + yjjd(pji) (3.12) 

-drf + {dx'^j + dyj + yfdipj) +2a(^^- 1^ drj dxj, (3.13) 



where we have dropped subleading terms suppressed by 1/a^ ~ exp(— 2A^) and 1/a ~ 
exp(-2iV* - 27V). 

Boundaries of the regions: The boundaries of regions 1, 11 and 111 arc defined by null 
surfaces originating from the surface {t = 0,x = 0,y,ip), propagating in the (r, a;)-plane. 
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This implies ds^ — with dy — d(f — 0, from which one can solve 

2 2 ^ ^ 

as the equation for the boundary between regions I and III; the approximation used was 
QrX ~ e^e^*~^ — e^* ^ 1. Since spacetime in region II is also RD-FRW with time 
coordinate th, the boundary between regions II and III is given by th ~ ^a^xjj. 



3.3.1 Region III 

We now define the coordinates and the metric in region III. We do this by first constructing 
a physically motivated ansatz for a coordinate system and interpolated metric. Then we 
vary this metric and interpolation, and compute certain coordinate invariant quantities to 
show that our ansatz produces the physically most reasonable interpolation in region III. 

We first fix the direction drjjj by using the symmetry of the spacetime by setting 

drjjj OC drj + drjj OC drj - l^^^i- (3-15) 

By orthogonality this fixes the coordinate xm to be of the form 

dxiii OC ^d-Ti + dxi, =^ xiii = h y^'Ti + xij , (3.16) 
for some function h. 

Since regions 1 and 11 are causally disconnected, identical RD-FRW spaces (with a relative 
boost between them), we expect evolution in each region to be identical. That is, we expect 
an event (e.g. recombination) that occurs in region I at r/ = tq will occur in region II 
at Tii — Tq. We are therefore motivated to take an ansatz for region III that matches 
from a point P — {tj — tq^xi — ^^^^) on the boundary between regions I and III to a 

point P = {tjj = To,xjj = ^^^^) on the boundary regions II and III along a surface of 
constant time {tui)- sec Fig. 2. Naively, inside region III this surface could take any shape, 
though clearly there are choices that are more physically motivated. The simplest physically 
motivated ansatz we can take is to demand the surfaces of equal time are linear in (t/,^/) 
(and (t//, xii)) coordinates. We will first make this choice, and then work out a more general 
ansatz in Sec. 3.3.2. We will show that to the order we are working in, all of these choices 
give similar results. 

The slope from P to P is independent of tq and this fixes 

Tin^ firi+^x^, (3.17) 

for some function /. Finally, we can fix the functions / and h by choosing (without loss of 
generahty) that the coordinates are continuous across the border so that points P and P are 
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given in Ill-coordinates as P — (tjjj — Tq.xju — :^^^) and P — [tuj — Tq.xih — — 9 
With these constraints the inverses of the functions are fixed as 

r\s) = «+^V^> (3-18) 

h~^{s) = s + sign(s)-s^ (3.19) 

where s is a dummy variable, and the sign(s) ensures that h is an odd function in s. These 
can be inverted to find / and h, but for our purposes these inverses are sufficient. 

Note that the coordinate system thus defined respects the symmetry of the spacetime; 
in Il-coordinates we have 

Tin = f(Tn + ^xny (3.20) 

xiii = -h (^■^'^ii + ^11^ ■ (3-21) 

Using the coordinate relations above and the RD-FRW metric (2.8), we can find the 
metric in regions I and II in Ill-coordinates. We compute these metrics at the boundaries'^ 
of region III, and then define the metric in region III as a linear interpolation in xm on 
slices of constant Tm between the boundaries. 

Using the relations (3.16) and (3.17), and the metric (2.8), we can write the metric at 
point P as 

f A \ Aa} 
dsjn\p = -drill + IH ) dxjjj + a?{dy]ii + Vnid^ni) ^dTnidxni, (3.22) 

where here a = \/Ctiii, and we have again neglected subdominant terms suppressed by e~^, 
e~-^* and (Ae)'^. Similarly, we can find the metric at the opposite point P 



2 



dsjiilp = -drjjj + a I IH v2to j dxjjj + a {dyjjj + ynid^ni) H ^driiidxiu. (3.23) 

Since P and P were arbitrary points on the boundaries, these relations specify the metric 
everywhere on the boundaries of region III. 



The interpolated metric: Wc define the metric in region III as a linear interpolation in 
xiu on the constant timeslicc Tm = Tq between (3.22) and (3.23). Since gi/J{P) = gll^{P) 
the linear interpolation on the timeslice Tm — Tq between P and P gives 

gii\To,xjn)=g'4'(P) = ~l. (3.24) 

^Note that going from region II to region III, the x-coordinate changes sign; Xn —Xm. This is because 
xii increases towards the left, and xm increases to the right. 

^°By 'metric on the boundary' we simply mean the values of the 4D metric tensor g^v on the boundary, 
not the induced 3D metric on the boundary. 

^^While we write g^-l^ for brevity, we of course mean g^^j^^^^^, and likewise for the other components. 
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Likewise, gJJ ^ Qyy and g^^ have the same value at P and P, and hence we fix those com- 
ponents to be independent of xm on the timeshce Tm — Tq. 

The remaining component is g^^l! , which changes sign between points P and P. The 
hnear interpolation then defines g^^J on the timeslice Tm — Tq &s 

gZ\ro,xni) ^ gii\P) ■ (^7=^^) = (3-25) 

which reduces to the correct boundary values as xjjj — )■ ±-^^^ = 

We have completely defined the interpolated metric on the timeslice Tm = Tq, and since 
To was arbitrary, this defines the metric in all of region III as 

dsjn = -dTijj + a^ ^1 + ^y/2Tni^ dx]ii + a^{dy]ii + y]jjdip]jj)- AaxnidTnidxni- (3.26) 



3.3.2 A general coordinate system 



The previous ansatz was the simplest physically motivated one. However, one could argue 
that a more general interpolation could capture the physics in region III better. In this section 
we will find a more general coordinate system, constrain it by physical considerations and 
then show that the effect on the galaxy fiow is minimal between these choices. We vary our 
coordinate system by defining 

T{Tin,Xni) = rjii + A6i{riij,xjjj), (327) 
X{Tni,xiii) = xni + A62{Tiii,xni), 

where di and S2 are a priori arbitrary, but will be constrained by physical considerations. The 
factors of A have been added for later convenience. Without loss of generality we pick 61 and 
62 to vanish on the boundaries of region III, to keep the coordinates continuous across the 
boundary. With this choice, it can be shown that the magnitudes of the Ss are e and 

|52| < e^. If (5i and 82 were bigger, T and X would not be timelike and spacelike respectively 
for generic Si and 62.^^ 

Using (3.27) and (3.16,3.17) wc solve for tj{T,X) and xj{T,X), and write the metrics 
in regions I and II in (T, X)-coordinatcs. As above, we then interpolate linearly in X across 
region III, giving the metric in region III as 

dsjjj = -dT^+a^{T) (1 + —V2T{1 - 6^)] dX^W{dyjjj+y'jjjdipjjj)-AaX{l-S^)dTdX, 



with a(T) = Vcf, and we have defined 

= -sign{X)2ajdTS2{T,X), 



(3.28) 



L ^ V^dxS2{T,X). (3.29) 
^^In the above we've assumed that Si and 62 depend polynomially on the variables, so that O(dTSi) ~ 



0(5^) etc. While more general forms could be considered, they do not seem physically well motivated. 
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The appearance of sign(X) above can be traced back to the the sign change xu — > —xm 
across the boundary. Due to the considerations above we generically have 
Note that 6i does not appear in the metric to this order, and thus we are free to set it to 
vanish, giving T = tjjj. 

We must demand that the boundaries T = ^a^X^ between regions I and III (and II 
and III) are null in this metric. The boundary is generated by the vector o^XOt + dx, and 
requiring this vector to be null gives 

S^{T,X)\,^S^{T,X)\,, (3.30) 

where the b indicates that this relation is valid on the boundary. 



Curvature in region III: At this stage, we have some constraints on the 5s, but a lot 
of freedom remains. We can ask what physical expectations we have about the behavior 

of physical quantities in region III and use that to further constrain the functions. For 
example, if the 5s were wildly varying functions, then the curvature in region III would be 
large in places, varying on small scales and possibly even negative. This behavior is clearly 
unphysical. We are thus motivated to examine the curvature scalar along a constant time 
surface to aid in constraining the 5s. We will find that with this criteria as well as symmetry 
considerations there are two natural choices, and the difference in fiow velocities between 
the two choices is at most 0(20% — 30%). 

The curvature scalar in region III is 
R = 3 (3 + 3aldT52 + ttr (3arX - A\/Yf^ 8x8x52 + 2Tal8^52 + 2Ta^ (^a^X - \/2t) 8^8x52^ 



ar{2T)2 

(3.31) 

The constraints on 52 to a large degree fix its functional form. By continuity of the coordi- 
nates 52 must vanish as X ^ ^"^^ symmetry of the spacetime it must vanish 
as X — )■ 0. Furthermore, as we have ^2 ~ e^, a fairly general^'^ polynomial function we can 
write down that satisfies (3.30) and the boundary conditions is 



where F is an unknown function of at most order unity, and ~ ^(1)- 



At the boundary between regions I and III we have = 1, and hence we expand 



2T 



F{s) — Yli ^n(l — •§)"■ for some constants We are mostly interested in geodesies relatively 
close to the boundary, so to a good approximation we can truncate 

^ 60 + &i fl - ^ V (3.33) 



2TJ V V2T/ 



^^Morc generally, a function of the form (^^^ j — ^"^^ , with 2c + b — a = 2 and b odd would 

work. Our choice, a = b = c = 1, seems physically the best motivated. 
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Figure 3: Graph of the value of the Ricci scalar for the two proposed interpolated metrics. 



In regions I and II (|X| > -^^) the spacetime is RD-FRW and the Ricci scalar vanishes. For 
the (t///, X7//)-metric the scalar is a positive constant proportional to e in region III given 
by the solid blue line, and for the (6o = — &i = ~io) system discussed in the text the Ricci 
scalar is given by the dashed purple curve. 



While this is clearly not the most general form of F in the interior of region III, it captures 
the physics close to the boundary, and by a suitable choice of 6o and hi we ensure that 
physics in the interior is also reasonable. We fix the constants 6o and hi by demanding that 
the curvature scalar vanishes on the boundary of the regions, and that the curvature scalar 
is differentiable at X = 0, as there are no sources along X = 0. Plugging (3.32) into (3.31) 
we find that these constraints imply 

&o = -fei = -^. (3.34) 

Note that the choice 6o = = gives the (r///, xju) coordinate system and corresponds to a 
constant curvature across region III proportional to e. In Fig. 3 we plot the scalar curvature 
on a constant time slice across the spacetime, for both h^ = hi = Q and 6o = —hi = —jq- 

In the next section we will show that the matter flow is largely unaffected between these 
two cases; the flow velocities in the 6o = —hi = are roughly 20%-30% higher than in the 
6o = &i = case. Thus our conclusions are fairly robust under the choice of the interpolated 
metric. 



4 Galaxy motion 

In the previous section we found a physically motivated metric for region III. Together 
with the metrics in regions I and II and the various coordinate relations we have a global 
description of the spacetime background to order e. In this section, we will find the motion 
of galaxies through this spacetime and show that the collision generically leads to peculiar 
velocities in the affected region. 
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In our approximation we will assume the matter that makes up galaxies was created 
at reheating, and was initially moving orthogonal to the reheating surface, i.e. initially 
generated by drj^ and drjj in regions I and II respectively. We track the future evolution of 
these geodesies and analyze the motion of the galaxies as seen by present day observers. 



4.1 Timelike geodesies in region III 

We are interested in observers originating in region I, who have entered region III in relatively 
recent times^'^. We define the time of the observer to be Tm = Tn ^ 0{e^^*). This implies 
the observer is at Xni = ~ 0{e^*~^); for larger values of x the observers will not yet see 
the collision. We will find the timelike geodesies in region III for this range of parameters. 

Galaxies that originate in region I follow the geodesies with constant Xi, yj and (pr while 
in region I. Upon entering region III they have four-velocity 

uo = drj ^ uldrjjj +u^d^jjj, with (4.1) 
ul = l-0{e'), < = ^(l + 0(e^)), (4.2) 
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where Uq = u^"' and = Uq"' , and Wq, ~ 0. 

We are thus interested in timelike geodesies in region III with initial velocity given above. 
Since the velocities are non-relativistic^^, the normalization condition \u\^ — —1 gives 



l-0{e'), ^T = s{l + 0{e')), (4.3) 



where s is the affine parameter on the geodesic. Making use of (4.3), we find the geodesic 
equation for to first order in e to be 

d'^xiH ^ 1 dxin Axni _ ^ 
drfjj Tin drill 2ar{2Tiii)i 



^^That is, we are most interested in observers who see a hot or cold spot that approximately takes up at 
most half of the CMB sky. A collision that takes up more than half of the CMB sky can be mapped in our 
approximation to one taking up less than half the sky by a boost. The former "spot" maps to the CMB rest 
frame and vice versa. To remain consistent within our approximations we will require that the CMB spot is 
less than about 80°, otherwise corrections of order can become important. It was argued in [17] and [34] 
that collisions that take up the entire sky are not observable, at least in the CMB temperature map as they 
appear as a pure dipole only. 

15A 

more careful computation yields ~ u^e^ <C tt^, and hence motion in the y-direction is negligible in 
the order in which we work. 

one does a similar computation for the CMB photons, they generically already have peculiar physical 
velocities of order 1 to satisfy the null condition. Therefore, the corrections due to the backreaction are 
subleading in e. If one computes their effect on the CMB temperature map as in [17] the effect from the 
backreaction multiplies a quantity that is already subleading, thus causing the effect to be 0{e^). It is 
therefore consistent at this order to ignore the backreaction effect for the CMB photons. 

^''The velocities given are in comoving coordinates; the corresponding physical three-velocity is w « 
u^a(r) <C 1 for relevant times r, hence non-relativistic. 
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This can be solved in the relevant range of coordinates to give 



X 



geod. 
Ill 



(tiii) = Xo + 



Am 



Tin 

To 



~2A 



In 



Till 
To 



+ 0{e'), (4.5) 



where {tq,Xq) is the point where the geodesic intersects the boundary between regions I and 
III, as shown in figure 5, and Uq is the a;-component of the four-velocity at that point. 



Peculiar velocity: Using the initial velocity found in (4.1) and the solution for the geodesic 
(4.5), we find that at time tjjj, the velocity of a galaxy which entered region 111 at time Tq 
is given by 

A 



dr. 



Ill 



2al 



To 

Tin 



(4.6) 



Note that by velocity we mean velocity in coordinate system III. To convert these to physical 
velocities, one must take into account the velocity due to the expansion of the universe, i.e. 
Hubble's law and convert to the CMB rest frame. 



The general coordinate system: The procedure above can be repeated for the (T, X)- 
coordinate system. In this case, the geodesic equation becomes 



d'^X 1 dX 
+ 



AX 



dT^ TdT 2a,(2r)i 



l + 4(6o + 6i) 



arX 



- 126i 



ttrX Y 



2TJ 



0, 



(4.7) 



and the initial velocity when the galaxy enters region III in these coordinates becomes 

A 
2a 



,x 



(4.8) 



The geodesic equation can be solved perturbatively to first non-trivial order in e to give 



AT. 




hi 



Tn 



T 



T 



4, 



To 
T 



(4.9) 

which reduces to (4.5) as 6o;^i 0. Here {Xq^Tq) is the point where the galaxy enters 
region III. Differentiating, we find the velocity along the geodesic as 



^jS^geod. ^ 



dT 



2al 




2{bo + h 



2^4 



(4.10) 



ispor galaxies that originated in region II, the analysis is identical aside from a straightforward sign flip of 
the second term. A corresponding change is made for subsequent formulas. It turns out that our results for 
the observable quantities we compute later on aiitomatically take this into account and so we suppress the 
further analysis of the region II galaxies here to simplify the presentation. The physical results we compute 
later are completely general. 
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Figure 4: Effect of the choice of coordinates on the flow velocities in appropriately normalized 
units, in the rest frame of the CMB (see text for discussion of the rest frame). The {60 = 
bi = 0} solution is portrayed with a thick blue line, while the {60 = —bi = —jq} case is 
shown with a dashed purple line. For t < tq the galaxies are in region I and hence have no 
flow velocity. 



To assess the effect of the choice of the interpolating coordinate system on the matter 
flow, we convert the velocity above to the (approximate) CMB rest frame^^ which is given 
by the coordinates of region I using the coordinate transformations (3.27), (3.17) and (3.16). 
This gives 



A 



2a2 



1 - 




2(60 + bi 




(4.11) 



= case. 



and in Fig. 4 we plot these velocities for the two preferred sets of values, bo = bi = and 
60 = —bi = —jq- Both sets lead to comparable flow velocities, with the 60 = —bi - 
generally leading to velocities that are roughly 20-30% higher than the bo = 61 
This shows that our model is fairly robust under the choice of a (physical) interpolating 
metric. For this reason, and for computational simplicity, we focus on the 60 = ^1 = case 
from now on; the results obtained for the other choice of interpolation will be similar, but 
would require numerical analysis. 



^^We are focusing on collision scenarios were the affected CMB region takes up less than half the sky 
(again if it takes up more, we can map it to a scenario where it takes up less than half by a boost), i.e. the 
observer's CMB sky is dominated by photons that originate in region I. Because of this, the CMB rest frame 
is that of region I to a good approximation. 
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4.2 Past lightcone of observers in region III 



In the previous section we found a general expression for the pecuhar velocity of galaxies 
relative to the CMB rest frame (4.11). We would like to find the observable consequences of 
these peculiar velocities, i.e. how big a patch of the sky they take up at a given redshift and 
how large the velocities are for reasonable scenarios. In this section we will work out how 
the size of the disk on the sky that contains galaxies with non-vanishing peculiar velocity 
varies as a function of A^, A^*, the redshift to the galaxies in question and also in terms of the 
effect on the CMB sky from the collision found in [17]. To do this, we take an observer who 
has recently entered region III, project back his past lightcone to a given redshift and then 
find the galaxies behavior within this lightcone. We will focus on redshifts of approximately 
0<^<1. 

Consider observers at point (r^, Vn, fn) in region III (and region III coordinate values), 
travelling along one of the geodesies found above (4.5). We wish to analyze the motion of 
galaxies as seen by these observers. As neither the metric nor the scalar field depend on y up 
to order C(e^), we are free to choose the origin of the {y, ip) plane to coincide with the location 
of the observers, so that y„ = 0. To find the past light cone, we define a coordinate system 
centered on the observers, (robs, Cobs, Uohs, 'Pobs), related to the usual region III coordinates by 

Tohs = Tin, Xohs = Xiii - Xn, (4-12) 

Vohs = yiii, (Pohs = (Piii- (4.13) 

We can find the null geodesies projected back from the observer in a similar manner to 
the procedure for finding the timelike geodesies of the previous sections 

Tobs(eobs) = I {V2^ - arU^Y (1 + ^iM) , (4-14) 

where ^obs = V^lhs + vlhs- Note that the leading term gives the parabolic past hghtcone of 
RD-FRW space. Using the leading term of the past lightcone is sufficient for our purposes, as 
the velocity (4.6) is already linear in e, and thus any subleading corrections to the lightcone 
will lead to sub-subleading corrections to the velocities of the galaxies, and will therefore be 
negligible to our purposes. The lightcone is illustrated in Fig. 5. 

The lightcone can be written in a more convenient form. Consider observers looking 
in a fixed direction 6'obs in the night sky, independent of (fobs- This translates to the III 
coordinates as tan^obs = = xni-x ~ constant, which implies ym — tan^obs(3^m ~ ^n)- 
Using this, the lightcone (4.14) can be written as 

XlII -Xn= \ COS 6'obs- (4.15) 



The collision disk: In [17] it was shown that a bubble collision generically produces a hot 
or cold spot on the sky. It is useful to ask how the galaxy flow will depend on the size and 
temperature profile of this spot. To that end, we express our results in terms of the angular 
radius of the affected disk on the CMB sky, given by tt — Oqub-"^^ The intersection of the 

^"Recall that in our coordinates 9 = it points towards the colhsion 
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two sloped parts of the reheating surface is given by Tm = xui = 0. Plugging this into the 
lightcone (4.15) gives the size of the collision disk seen by the observers as (see Fig. 5) 



cos6'cMB = 




(4.16) 



Observable galaxies: Consider other galaxies as seen by our observers; illustrated in 
Fig. 5. The intersection of the geodesic of the other galaxy with the past lightcone of our 
observers is given by the equations 



Galaxy geodesic: 



Xg Xgg 



Lightcone: x„ — 



At. 



90 




1 



COS 9, 



obs; 



(4.17) 
(4.18) 



where {Tg,Xg) is the point of intersection and (T^ojXgo) is the point where the other galaxy 
crossed the boundary between regions I and III; refer again to Fig. 5. 

The peculiar velocity of the galaxy (in region III coordinates) at the intersection is given 
by (4.6) Vgii — V^go/^s- To write this in variables more appropriate to the observer, 
we first use the equation for the geodesic of the galaxy to solve for Tg^ = Tg^^{Tg,Xg). As 



(r^Q, Xg^^Ug^) is on the boundary between the regions, we have Xg^^ 

^2 ''go AXn 



'go 



-, and we get 



Xq -\- 



(4.19) 
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We insert this into the velocity, and ehminate Xg in favour of ^obs and t„ using the 
equation for the hghtcone. Furthermore, as hght in FRW space redshifts as oc 1/a, the 
redshift z is given by 1 + 2; = The velocity to leading order is 

= ^ [- cos ^cmb(1 + z) + z cos ^obs] , (4.20) 

as a function of the direction in the night sky that is being observed, the angular radius of 
the CMB coUision disk, and the redshift to the observed galaxy 



A 

2^ 



1 + z) -\- zcosOobs 



The CMB frame: The velocity computed above is in Ill-coordinates; it can be trans- 
formed back to I-coordinates, making up the (approximate) rest frame of the CMB, using 
the coordinate relations (3.17) and (3.16). This gives the velocity relative to the CMB rest 
frame as 



^g'af = ^gal' - ^ + ^(^)) - [1 + COS^cMb(1 + z) - Z COS ^obs] • (4.21) 



r 



The relation above is valid in region III; galaxies in region I are moving along drj and thus 
have no velocity with respect to the CMB frame. We can use this to find the size of the disk 
of affected galaxies by setting the velocity Vglf to vanish for a fixed redshift z, in which case 
(4.21) implies 

cos6lgai(2;) = - (1 cos6Icmb(1 + z)) = cos6'cmb + ^^^^ . (4.22) 

z z 

All the galaxies inside this disk, ^^obs > ^gai, are inside region 111 and thus are affected by the 
collision. Note that 6'gai < 6'cmb^^, and hence the disk of affected galaxies is larger than the 
collision disk, which is also clear from Fig. 5. Inside the disk of affected galaxies, we can see 
that the velocity at a fixed redshift decreases linearly with cos ^obs from its maximum value 
to zero at the edge of the disk. 



5 Results 

In this section we examine the effect of a bubble coUision on the large scale galaxy flow in 
two scenarios: 

1. A coUision that produces a small spot on the CMB sky, akin to the WMAP cold spot. 

2. A collision that produces a larger disk on the sky. 

We will find that making our results agree with CMB observations will constrain our possible 
flows. 



^^This formula is general and therefore also valid for galaxies that originate in region II, taking into account 

the sign flip pointed out in footnote 18. 



^^Recall again, the angular radius of each disk is given hy n — 6 respectively. 
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Physical velocities: In the previous section we found the comoving pecuhar velocities of 
galaxies inside the observer's past lightcone (4.21). To convert them to a physical velocity 
that can be observed we multiply by the scale factor at the time the galaxy intersected with 
the observer's past lightcone to get (we work in units where the speed of light, c = 1, and 
our expression is vahd to 0(e)) 



phys 



, , A a(r„) / 1 + cos 6'cmb + ^(cos 6'cmb - cos 6'obs) \ 

^ = ) 

1 + COS 6'cmb + ^(cos 6'cmb - cos 6'obs) \ . n,-n /r i \ 

where we wrote a{Tg) in terms of the redshift and a(T„), and in the last line we used ^izii) i — 



Effect on the CMB: It will be natural to express the galaxy flows we find in terms of 
the corresponding effect from the collision on the CMB itself. This effect was found in [17] 
and we refer the reader there for details. The redshift seen by an observer is defined as the 
ratio of the energies of the photon as it was emitted and absorbed, 

-E'n Cl^^r) . . 

Er a{Tn) Tr 

For the unperturbed part of the reheating surface this is simply given by = ar/an = 
g-v, _ 7"0/x^. In the direction of the collision, 6^obs > 6'cmb, reheating happens slightly later 
(earlier) if ^4 > (^4 < 0) (see Fig. 2) and the redshift is given by 

e-^Vl-2Ax,f« e"^* - [cos^obs - cos^cmb]^ 

(5.3) 

Here < is the a:-coordinate of the intersection of the observer's past lightcone with the 
perturbed reheating surface, in a given direction ^obs, and we have expressed Xj. in terms of 
more physical quantities using the equation for the lightcone, (4.15). The variation from the 
average temperature due to the collision is^^ 

— = (cos 6'cmb - cos 6'obs) = (cos 6'cmb - cos 6'obs) (2^e^*"^) , (5.4) 

and as AT/T° is restricted by experiments, this expression can be used to give the matter 
flow velocities in terms of experimental data. 




^"^There will of course be the usual Gaussian fluctuations in the temperature map. They are not important 
for our analysis here, but a full treatment taking them into account is given in [17]. 
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Figure 6: Flow velocity for the small collision disk of scenario 1 at redshift z = .02 as a 
function of the angle on the sky. The velocity profile for general redshifts looks similar, 
always decreasing from the same maximum value to zero linearly in cos 6'obs to the edge of 
the disk. Outside the disk the peculiar velocity always vanishes. 



The largest flows and temperature variations are always in the direction of the collision, 
5'obs = tt), in which case we get 



Ae 



N,-N 



AT 



n / max 



2 1 + COS 6'cMB 

We can use this to express the flow velocity in terms of the temperature of the spot 



(5.5) 



1 /AT 



^phys 



2 V TO 



n / max 



1+Z 



l + Z 



COS 6. 



CMB 



cos6', 



obs 



1 + COS 6'cMB 



(5.6) 



The largest flow we can measure is 
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AT 



2^0 



(5.7) 



note that this is independent of the size of the collision disk and the redshift. 



5.1 Scenario 1 : A small collision disk 

We first examine the case of a small collision disk. If we take the WMAP cold spot as a 
rough guide, it is ~ 73fiK colder on about 5° scales [21, 22, 23, 24, 25]. It was shown in [17] 
that the effect on the CMB of a collision is generically a spot, with a temperature profile 
that is linearly decreasing in cos^obs from the center of spot to its edge, hence its average 

^^Of course, there will be large fluctuations based on random local behavior for a given galaxy cluster. 
What we mean here is the maximum coherent contribution from a collision in our scenario. 
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Figure 7: Angular radius of the disk of affected galaxies for the small collision disk example 
of scenario 1 . The solid blue line is the size of the galaxy disk and the dashed purple line gives 
the (constant) size of the cold spot for comparison. Note the galaxy disk can be significantly 
larger than that of the CMB spot, especially at low redshift. 

temperature is ~ Tspot/2, where Tgpot is the temperature at the center of the spot. To give 
us an idea of the type of flows a collision that causes a spot like this generates, we take 
T° — Tspot = 300 fiK and angular radius 180° — ^cmb ~ 7°. This gives a flow velocity profile 
given by (5.6) with {AT/T^)^^^ ^ 5.5 x 10~^. The maximum detectable flow velocity is 

\v^^^,\ ~ 0{5 ■ 10-') ~ 0(15km/s). (5.8) 

At a given redshift z, the magnitude of the velocity inside the disk of affected galaxies 
will decrease from the redshift independent maximum value jf^y'^l to zero linearly in cos^obs 
at the edge of the disk, we display a typical profile in Fig. 6. The size of the disk is given 
in (4.22). In Fig. 7 we plot the angular radius of this disk as a function of redshift. We see 
that the angular size of the galaxy disk is always larger than that of the corresponding CMB 
spot, particularly at lower redshifts. 



5.2 Scenario 2 : A large collision disk 

Here we examine a larger collision disk, but with subsequently smaller AT so that it can still 
agree with present observational constraints. Such a disk could still contribute in the CMB 
data, for such effects as the hemispherical power asymmetry [26] or spatial correlation of the 
lower multipoles. We take cos6'cmb = —0.3, corresponding to a disk of angular radius ~ 70° 
on the sky and a corresponding (AT/T°)max = 5 x 10~^. We then find a disk of affected 
galaxies of angular size in Fig. 8 and maximum flow velocity 

|r;Xsl~(^(10km/s). (5.9) 
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Figure 8: Angular radius of the disk of affected galaxies for the large collision disk example 
of scenario 2. The solid blue line is the size of the galaxy disk and the dashed purple line 
gives the (constant) size of the CMB spot for comparison. Note the galaxy disk can be 
significantly larger than that of the CMB spot. At lower redshifts we observe a flow across 
the entire sky. 

In Fig. 8 we see that for redshifts out to 2; ^ 0.6 we observe a coherent flow across the 
entire sky. Even out to a redshift of one we still see a significant fraction of the sky flowing. 

6 Conclusions and comparison with observation 

In this paper we have shown that bubble collisions generically lead to a large scale bulk flow of 
galaxies. There is a disk on the sky of affected galaxies, always larger than the corresponding 
spot for the CMB temperature map, that decreases in size with redshift (though it may take 
up the entire sky for all redshifts observed). Within this disk, galaxies will coherently flow, 
with a peculiar velocity that decreases in magnitude from the center of the disk and is 
proportional to the corresponding temperature at the center of the spot in the CMB map. 
This is a distinctive signal profile and is in principle measurable with improved observations 
or analysis of current cluster catalogs. It would be very interesting if a bulk galaxy flow 
could be correlated with a corresponding disk or spot in the CMB data as this would be 
evidence for bubble collisions and the string landscape. 

Recently, two groups have claimed to observe bulk flows of galaxy clusters [20, 18, 19]. 
They both measure a flow of (9(400 — 800 km/s) out to a redshift z ~ 0.2. Both measurements 
are subject to a variety of random and systematic uncertainties (e.g. filtering methods and 
the conversion from temperature to velocity). They seem to have reasonable directional 
agreement, but the flow measured by Watkins et. al is considerably smaller than that of 
Kashlinsky et al. and they measure over different redshift ranges (Watkins et. al catalog is 
dominated by objects with z < 0.03 and Kashlinsky by objects with 0.04 < z < 0.2). The 
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flows we find appear to be an order of magnitude smaller, as our fiows are constrained by the 
corresponding effect in the CMB. Larger flows would lead to a very cold or hot spot in the 
CMB, likely ruled out by current observations. In addition, since our peculiar velocity profile 
falls off the further we are observing from the center of the disk of affected galaxies, and 
the disk itself may not take up the whole sky, the corresponding average peculiar velocity 
measured would be smaller. However, we have made several approximations which could 
affect the measured fiow: 

• Wc have ignored any self-interaction and gravitation of the radiation and matter fiuids. 
It is possible that as the fluids clump, this could affect the backreacted metric in region 
III and hence the flow. 

• We have ignored the transition to matter domination. Taking this transition into ac- 
count would likely increase the flow for a subtle reason. In a matter dominated universe 
peculiar velocities will redshift more quickly. This implies that galaxies entering region 
III will transition more quickly to flowing along drjjj- Transforming back to the CMB 
rest frame this results in higher measured peculiar velocity with respect to the rest 
frame. 

• The small freedom we have in choosing the metric interpolation in region III can have 
an 0{1) effect on the ffow. The flows computed above were for the choice that seems 
to produce the smallest flow. Other choices increase the flow velocity. 

• We have likely not accounted for many more detailed effects of the cosmological evo- 
lution that could affect the flow. 

Because of these approximations, it is not hard to envision a more detailed treatment flnding 
an enhanced flow. 

Many open questions on bubble collisions remain. It would be interesting to carry out a 
more rigorous model of cosmological evolution and perturbation theory in this background 
and flnd a more accurate picture of these effects. In this paper and [15, 17] we have mapped 
out two possible effects of bubble collisions resulting from the string landscape, namely the 
effect on the CMB temperature map and on galaxy cluster motion. One could search for 
other effects, e.g. polarization of the CMB as our toy models bear some relation to those of 
[35] and it may be possible to further correlate these signals. This effect could be particularly 
interesting in light of new polarization data expected from the Planck experiment. One could 
also look for signals in 21cm radiation data. On the observational side, most of the ciurcnt 
measurements of the CMB have been presented in momentum (£) space, yet these effects are 
most apparent in real space. Each collision would lead to a disk on the sky and one could 
search the CMB and other data for traces of disks, some of which might overlap. While 
many of these things are likely to be difficult or require improved measurements (e.g. 21cm 
telescopes), we feel that cosmology and astrophysics offers a unique opportunity to possibly 
detect effects from high energy physics and string theory. 
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